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ANALYSIS OF TRANSONIC FLOW ABOUT LIFTING 
WING- BODY CONFIGURATIONS 

Richard W. Barnwell 
Langley Research Center 

SUMMARY 

An analytical solution is obtained for the perturbation velocity potential for tran- 
sonic flow about lifting wing-body configurations with order-one span-length ratios and 
small reduced- span-length ratios and equivalent- thickness-length ratios. The analysis 
is performed with the method of matched asymptotic expansions. The angles of attack 
which are considered are small but are large enough to insure that the effects of lift in 
the region far from the configuration are either dominant or comparable with the effects 
of thickness. The modification to the equivalence rule which accounts for these lift effects 
is determined. An analysis of transonic flow about lifting wings with large aspect ratios 
is also presented. 


INTRODUCTION 

It is well known that the outer part of transonic flow fields about lifting configura- 
tions at small angles of attack is governed largely by the effects of thickness. It has been 
shown by Oswatitsch and Keune (ref. 1) that the flow in the region far from nonlifting wing- 
body configurations is mathematically equivalent to that about axisymmetric bodies with 
the same cross-sectional area distributions. Whitcomb (ref. 2) has demonstrated experi- 
mentally that the wave drags of wing-body combinations at zero angle of attack are the 
same as those of the equivalent axisymmetric bodies. The analysis of Heaslet and 
Spreiter (ref. 3) shows that these zero- angle- of- attack results also apply to configura- 
tions at angles of attack of the order of the wing thickness or less. The main purpose of 
this report is to present an analytical solution for transonic flow about lifting configura- 
tions at angles of attack large enough to insure that the effects of lift either dominate or 
are comparable with the effects of thickness in the outer region. 


The present problem and two similar problems which have been analyzed previously 
are compared in figure 1. The problem of transonic flow past slender bodies at angles of 
attack of the order of the equivalent-thickness-length ratio has been treated by Messiter 
= (ref. 4), Hayes (ref. 5), and Lifshits (ref. 6) with the method of asymptotic expansions. 

IT The problem of transonic flow past wing-body combinations with spans of order one at 



Previous studies Present study 



a o c Body thickness ratio a oc Wing thickness ratio a oc Body thickness ratio 
Figure 1.- Comparison of past and present transonic studies. 

angles of attack proportional to the wing thickness-length ratio or smaller has been stud- 
ied by Hayes (ref. 5) and by Cheng and Hafez (ref. 7) with the method of matched asymp- 
totic expansions and by Heaslet and Spreiter (ref. 3) with the integral equation formula- 
tion. Since the flow is thickness dominated for both of these problems, the equivalence 
rule of Oswatitsch and Keune applies. As is indicated in figure 1, the present problem 
involves flow about wing-body combinations at angles of attack which are small but much 
larger than the wing thickness. In the first attempt at solving the present problem, Cheng 
and Hafez (ref. 7) obtained results which indicated that when the effects of lift and thick- 
ness were comparable in the outer region, the expressions for each of these effects are the 
same as that obtained when the effects of thickness are dominant. In the preliminary ver- 
sion of the present treatment (ref. 8), it was shown that, to the contrary, when the effects 
of lift and thickness are comparable in the outer region, there is a source flow due to lift 
which is of the same order of magnitude as the source flow due to thickness. Later in 
this paper, it is argued on physical grounds that this type of influence by lift is reasonable. 
As a result of the source flow due to lift, the equivalence rule of Oswatitsch and Keune is 
not applicable to the present problem. It should be noted that the radial length scale in 
the outer region obtained in reference 8 for the present problem was different from that 
obtained for the thickness-dominated problem. This difference in radial length scale 
resulted from the choice of relationship between the angle of attack and the equivalent 
thickness ratio which was used. In reference 9 Cheng and Hafez show that there is a 
second relationship between the angle of attack and the equivalent thickness ratio which 
leads to the same source flow due to lift obtained in reference 8 and which yields the same 
outer-region radial length scale found for thickness-dominated flows. This second length 
scale is physically more realistic and thus is the one used in this report. 

In this report the analytical solution for transonic flow about lifting configurations 
with span-length ratio of order one and small reduced span-length ratios ^the product of 

the span-length ratio and the factor l/|i - m£| when Moo I s the free- stream Mach num- 
ber^ is derived in detail with the method of matched asymptotic expansions. This solution 
is determined to within an arbitrary additive function of the length along the axis which 
cannot be determined with the present method. There is a brief discussion of how this 
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iunction can be determined by use of numerical techniques. Also a study is made of the 


pr: problem of lifting transonic flow about configurations with large aspect ratios. 


SYMBOLS 


aESL. _ A (x) ,B (x) ,C 



(x) arbitrary functions in unmatched outer expansion 

functions in twist and camber potential (see appendix B) 
constant of proportionality in equation (54) for angle of attack 
wing semispan 
pressure coefficient 

constant of proportionality of order one (see eqs. (3) and (95)) 
exponent, 0 g c s 1 

Cartesian coordinates locating vortex core relative to wing tip (see fig. 8) 
complete elliptic integral of second kind 
nondimensional thickness distribution of body 
nondimensional thickness distribution of equivalent body 
nondimensional dipole strength distribution 

arbitrary additive functions in outer expansion (see eq. (68)) 


g(x,y) twist and camber distribution of wing 

g n (x),g m>n (x) arbitrary additive functions in inner potentials <£ n and $ m)I1 

gg (x) arbitrary additive function in inner thickness potential <p^ 

H(x) secondary source- strength distribution in outer expansion given 

by equation (51) 
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h(x,y) 

K 

k 

L 

a 

Moo 

m(x,y) 


thickness distribution of wing 


_ 1 - M£ . 
S 1 


e ; also complete elliptic integral of first kind 


modulus of E and K given by equation (26) 
characteristic length in x-direction 

characteristic length in x-direction, configuration length unless otherwise 
noted 

free- stream Mach number 


= -i(m + - m_) 


m + (x,y),m_(x,y) functions in equation (29) for potential 4> ^ near configuration 

surface 


N 


nonnegative number 


q x ,q r ,q 0 body-oriented cylindrical polar velocity components (see eqs. (2)) 
r radial body- oriented coordinate 


r b (x) 

S(x) 

si 


t 

tw 


radius of body surface 

Mirels’ S- function (also Mangler's H-function) 

lowest order source strength of configuration in outer region 
(see eqs. (69), (70), and (71)) 

distance in cross- flow plane along wing and leading- edge vortex sheet from 
center of wing 

maximum thickness of equivalent body 
maximum thickness of wing 


U D 


free- stream speed 
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u,v,w 


v c» w 


c 


w 

X 


x,y,z 


Y 

y 2 (x),y 1 (x) 

Z 


z w 

a 

P 

r v 

y 

y(x,s) 

A(p 

5 

e n 

Vl 

0 


body- oriented cylindrical polar perturbation velocity components 

body- oriented Cartesian perturbation velocity components in 
y- and z- direction, respectively 

complex potential 

complex variable, y + iz 

body-oriented Cartesian coordinates (see fig. 3) 
complex function defined by equation (A2) 
leading and trailing edges of wing 
transformed complex variable defined by equation (A16), rj + i£ 
z- coordinate of wing surface 
angle of attack 

angle defined by equations (A9) or (A14) 
strength of vortex core 
ratio of specific heats 

vortex strength of wing or vortex- sheet segment 
potential jump across vortex- sheet segment 
equivalent body thickness ratio, t /ft 
nth gage function in outer expansion for cp 
coordinates in transform planes (see figs. 8 and 9) 

V ~ coordinates of leading and trailing edges of wing (see fig. 9) 
body- oriented polar angle (see fig. 3) 


A 

A. 

Mn 


$ 


n 


0m, n 

0n 


= b. 
L 

_ b 

SL 


nth gage function in intermediate expansion for <p 
stretching parameter for outer radial coordinate 
infinitesimal radius 

angle in transform plane depicted in figure 8 
nth potential in outer expansion for (p 

potential in inner expansion for cp with gage function log™^Ljsin n 
potential in inner expansion for cp with gage function sin 11 a 


a 


^la’^lb components of potential 4>i (see appendix B) 
02a>02b’02e’02d components of potential <p 2 (see appendix C) 


<P 


perturbation velocity potential 


(Pq ,(p§ j thickness potentials in inner expansion for cp 


'•n 


nth potential in intermediate expansion for cp 


0 component of thickness potential cp^ given by equation (30) 

w,a^,a >2 angles shown in figure 10 


two-dimensional Laplace operator in cross-flow plane (see eqs. (12) and (21)) 


Superscripts: 


complementary solution 


particular solution 
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Subscripts: 


max maximum value 

v vortex 

A bar over a symbol denotes outer variable (see eqs. (9)). A tilde over a symbol 
denotes inner variable (see eqs. (8)). A circumflex over a symbol denotes an intermediate 
variable (see eqs. (72)). Primes denote differentiation with respect to x, x, or x. An 
asterisk denotes a complex conjugate. 

EFFECT OF LIFT AT NEAR-SONIC SPEEDS 

It is well known that the effect of thickness is to deflect streamlines outward from 
the configuration. It can be shown that for near- sonic flow, lift can also have this effect 
in addition to the usual downwash effect. The manner in which lift causes the outward 
deflection of streamlines is depicted schematically in figure 2. The cross-sectional area 


Stream tube size 
near minimum in 
transonic free 
stream 




Velocity increases above wing; 
stream tube size increases 



(I 




Velocity decreases below wing; 
stream tube size increases 



Figure 2.- Outward displacement of streamlines by both lift and thickness effects. 

of stream tubes is minimum where the flow is sonic. The effect of lift is to increase the 
fluid speed in stream tubes above the wing and to decrease the speed in the tubes beneath 
the wing. For configurations traveling at Mach numbers near 1, both the increase and 
decrease in fluid speed are deviations from near- sonic flow. Thus for sufficiently large 
angles of attack, the cross sections of practically all the stream tubes about the body 
increase so that the streamlines are deflected outward more than they would be by thick- 
ness effects alone. The magnitude of the angle- of- attack range in which this phenomenon 
occurs will be established subsequently. It should be noted that the phenomenon does not 
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occur for completely subsonic or completely supersonic flow, where an increase in 
stream- tube size on one side of the wing is compensated by a decrease on the other. 


SLENDER CONFIGURATION PROBLEM 


Wings are described traditionally as being slender if the reduced span (the product 

of the span and the factor ]J\l - m£| ) is much smaller than the length. In this section an 
analytical solution is obtained for transonic flow past slender configurations with span- 
length ratios of order one which are at angles of attack large enough to insure that the 
effects of lift and thickness are comparable at large distances from the configurations. 
The method of matched asymptotic expansions is used. The solution is determined to 
within an arbitrary additive function of distance along the axis which cannot be obtained 
with the present method. There is a brief discussion of methods which can be used to 
determine the additive function. 


Problem Description 

In this subsection the problem is described and the governing equation and boundary 
conditions are established. The problem is to analyze transonic flow about lifting wing- 
body configurations with small equivalent- thickness ratios 6 and semispan-length 
ratios X of order one which are at angles of attack a so that the effects of lift in the 
outer region either dominate or are comparable with the effects of thickness. Let the 
free-stream Mach number, the configuration length and semispan, and the maximum thick- 
nesses of the wing and equivalent body be M^, b , t w , and t, respectively. The 
flow fields which are considered in this section are characterized by the parameters 


|Moo - 1 | « 1 

X = £ = 0(1) 

6 = i « 1 ) 

¥« i 


( 1 ) 


a « 1 j 


Body- oriented cylindrical polar and Cartesian coordinate systems are used in this 
report. These coordinate systems are shown in figure 3. For a vehicle at an angle of 
attack a, the total velocity components in the x-, r-, and 0-directions are 
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( 2 ) 


q x (x,r,0) = Uoo cos a + u(x,r,6>) 

q r (x,r,0) = sin o sin 0 + v(x,r ,0) > 

q e (x,r,0) = Uoo sin a cos 9 + w(x,r,0) ^ 

where U m is the free-stream velocity and u, v, and w are the perturbation velocity 
components in the x-, r-, and 0-direction, respectively. The Cartesian perturbation veloc- 
ity components in the y- and z- directions are designated as v c and w c , respectively. 



Figure 3.- Frame of reference. 

Consider a wing-body combination composed of a slender axisymmetric body and 
a thin wing which passes through or near the body axis. It is assumed that the cross- 
sectional areas of the wing and body are, in general, of the same order of magnitude. 
Since the cross-sectional area of the wing is proportional to both t^ and bt w , the wing 
thickness can be expressed as 
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( 3 ) 



The body radius r b and wing thickness z w , which are depicted in figure 4, are specified 
by the equations 

r b = t F(x) ( 4 ) 

and 

z w = f sin o' g(x,y) ± t w h(x,y) (5) 

respectively, where the plus sign in equation (5) is used for the top surface of wing and 
the minus sign is used for the bottom surface. The function F(x) is the nondimensional 
radius of the body, and the functions g(x,y) and h(x,y) give the twist and camber dis- 
tribution and the thickness distribution of the wing, respectively. As shown in figure 4, 




(a) Planform. (b) Cross section. 

Figure 4.- Vehicle geometry. 

the leading and trailing edges of the wing are specified by the functions y 2 (x) and y^x), 
respectively. It is assumed that the wing twist and camber are scaled by the quantity 
l sin a (that is, the order of magnitude of the twist and camber is assumed to be l sin a) 
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so that the effects of these quantities are of the same order of magnitude as those of angle 
of attack. The exact boundary conditions at the body and wing surfaces are 


v(^x,r b ,0) = -Uoo sin of sin 



cos a + u 



t 


dF(x) 

dx 


(6a) 


and 

w c (x,y,z w ) = -U M sin a + ju^ cos a + u(x,y,z w )J 


H sin a 


3g(x,y) 


ax 


± fc 


9h(x,y) 


W 


ax 


+ v c (x,y,z w ) 


sin a 


Qg(x,y) . . 

ay * *w 


9h(x,y) 

9y 


(6b) 


respectively. 

It is assumed that the flow fields under consideration are isentropic. Although this 
assumption is not strictly valid when shock waves are present, it can be shown (ref. 10) 
that the largest term which is affected by the nonisentropic condition is of the order 
(l - M^,) , which is very small for transonic flow. Hayes (ref. 5) and Cole and Messiter 
(ref. 11; also ref. 4) show that to the degree of approximation employed in this report, the 
isentropic assumption is valid. 

If the flow is isentropic, the perturbation velocity potential cp can be defined so 

that 


u 


a cp 
ax’ 


V 


dcp 

ar 


w 


_ i»£ 

r a 0 


From equation (5.5) of reference 4, it can be shown that the exact partial differential equa- 
tion governing cp is 


i - Mi 


\ a 2 jP 
ax 2 ' 


a 2 <? l W 

ar 2 r ar 


_1 

r 2 a e 2 


— = 


!sin 2 a( ||: Si n e + l^cos B) 


0) + sin 2 0 - + (^-S- sin 6 + - 9 ^ cos sin 6 

/ dx z \ar z r 9ray 


(fe Sine + r^f cos ®)c° se + 


cos G 




cos G 


cos O' 
U* 


1/ , 

|(y + 1) ^r^2 


(Equation continued on next page) 


11 



2 a£_a> 

dr dxdr 


dtp 


8 2 t p 


dQ dxd6 


+ (r 


-D«& + m + 4ij + ®i2zL - nte 

; 8x * r ar + 7? se 2^ U OT [ y ^8r 


- sin 0 




r 3 80 


d 2 (p 

dB 2 


COS 0 


+ 



2 a^£ 8£ a 2 <y> 
dx dr dxdr 


2_acpd(p d 2 (p _2_ _a^ _a<£ a 2 y /a<^\ 2 a 2 <i y - l /a<^ 2 f a^\ 2 
r 2 ax 80 8x80 + r 2 3r 90 a r 30 + y 8 x J 0 x 2 + 2 |Vax/ + \9r/ 


. l /M 2 /j>V , a 2 » 

r 2^0/ _\ax2 ar2 


I 

r ar 



( 7 ) 


Application of Method of Matched Asymptotic Expansions 

This method is applicable to problems for which the characteristic length scale of 
disturbances in a limited region of interest, called the inner region, differs from that in 
the surrounding, or outer, region when one or more control parameters are small. As a 
result of this difference in length scales, the magnitudes of various terms in the governing 
equations differ in the two regions so that the approximate forms of the governing equa- 
tions in the two regions are different. The formal procedure is to determine asymptotic 
expansions in the inner and outer regions and to match these expansions term by term 
according to the matching principle of Kaplun and Lagerstrom (ref. 12) in the zone where 
the two regions overlap. This matching of terms serves to determine coefficients and 
gage functions in the expansions. If the inner and outer regions do not overlap, as in the 
case for the problem treated in this report, the formal procedure includes the determina- 
tion of a third expansion in an intermediate region which overlaps the inner and outer 
regions. This intermediate expansion is then matched to both the inner and the outer 
expansions. A detailed discussion of the method of matched asymptotic expansions has 
been given by Van Dyke (ref. 13). 

The problem of transonic flow past lifting configurations with span-length ratios of 
order one is well suited for treatment with this method since the magnitude of the length 
scale in planes normal to the axis is the same as that of the span near the vehicle but 
is much larger than that of the span at large distances from the vehicle if the reduced 

semispan- length ratio \/|l - m£| b fl is small. In this report the length scales for 
the axial and radial coordinates in the inner region are chosen to be the configuration 
length L and the semispan b, respectively. Therefore , the independent variables in 
this region are x, r , and 9 where 
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(8) 




The length scales for the axial and radial coordinates in the outer region are fL and 
L/v, respectively, where 

v « 1 



Thus, the independent variables in the outer region x, r, and 9 are 







where \ = j. 


(9) 


It will be shown that as indicated previously, the regions of validity of the inner and 
outer expansions for the present problem do not overlap. Consequently, an intermediate 
expansion is obtained and matched with the inner and outer expansions in order to show 
that the formal matching procedure of Kaplun and Lagerstrom (ref. 12) can be applied to 
this problem. However, it is first shown that the inner and outer expansions can be con- 
structed with a less formal and also less complicated procedure. This procedure con- 
sists of the determination of the gage function of the leading term of the outer expansion 
in terms of the stretching parameter v, the determination of the inner expansion for 
large values of the inner radial variable r, the transformation of this expansion to outer 
variables, and the matching of the gage function of the leading term in the transformed 
expansion to the known expression for the gage function of the leading term in the outer 
expansion. 


Scaling of Basic Parameters 

The relationships between the parameter v, the free- stream Mach number M^, 
and the magnitude of the velocity potential in the outer region can be established from 
physical considerations. Let the expansion for the velocity potential in the outer region 
be written as 


(p(x,r,9) = ^ ^(xjr,#) 

i=l 


(10) 
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where 


ej « 1 


e i+l <<: e i 


From equations (7), (9), and (10), it can be shown that the first approximation to the gov- 
erning equation in the outer region is 


(l - + 1^1*1 = (r + 


(ID 


n 

where Vg, the Laplace operator in planes normal to the x-axis in terms of outer vari- 
ables, is written as 


y2 3 2 j 1 8 | 1 8 2 

2 ar 2 r 9r r 2 90 2 


( 12 ) 


and where the prime denotes differentiation with respect to x. If equation (11) is to be 
valid for completely subsonic and completely supersonic flows, the magnitudes of the 
terms on the left-hand side must be the same. If the sign of the coefficient of the deriva- 
tive in equation ( 11 ) is to change, as it must for transonic flow, the magnitude of the 

term on the right-hand side must be the same as that of the first term on the left. Con- 
sequently, the quantities 1 - m£, v 2 , and in equation ( 11 ) must be of the same order 

of magnitude. In this report this fact is expressed as 


V 2 = = 


M 2 


K 


(13) 


where K is an order- one constant. Note that equation (11) is, in general, of the mixed 
elliptic-hyperbolic type. 

As noted in 'Introduction," the outer regions of transonic flow fields about slender 
bodies at angles of attack of the order of the equivalent thickness ratio and about wing 
bodies with span-length ratios of order one at angles of attack of the order of the wing 
thickness are dominated by thickness effects. Oswatitsch and Keune (ref. 1) show that if 
thickness effects are dominant, the radial stretching parameter v and the equivalent 
thickness ratio 5 are of the same order of magnitude. It should be noted that it is com- 
mon practice to equate these two quantities. In reference 9 Cheng and Hafez point out that 
it is physically realistic to assume that the relationship between u and 6 which is valid 
for thickness-dominated flows should also hold for flows where the effects of lift and thick- 
ness are comparable. Consequently, it is assumed in this report that v and 6 are 
related by the equation 

14 


v = 6 


( 14 ) 


if the effects of lift and thickness are of the same order of magnitude in the outer region 
and by the inequality 

v » 6 (15) 

if the effects of lift are dominant. 

It should be noted that it remains to establish the relationship between the angle of 
attack a and the other basic parameters. 


Inner Expansion 

The perturbation velocity potential in the inner region may be written as 


<P(x, r,0) = <6 2 log e ^^ 6>1 (x,r,0) + 5 2 <p 6 (x,r,0) + X sin a 0 1 (x,r,0) 


00 m 

+ ^ \3m-2 sin m a 0 m (x } r,0) + ^ iogg^jj) 0 m , n (x,r,0) 

m= 2 n= 1 


(16) 


The first two terms are the usual lowest order thickness terms, the third term is the 
usual lowest order lift term, and the remaining terms are the higher order lift, thickness, 
and lift- thickness interaction terms. The thickness potentials have the leading sub- 
script 6 , and the potentials with gage functions proportional to sin m a have the lead- 
ing subscript m. The potentials with gage functions proportional to l°Se(\^) h ave a 
second subscript n. Since a unique relationship can and will be established between the 
angle of attack a and the equivalent thickness ratio 6, it is possible to express all the 
gage functions in equation (16) in terms of only one of these parameters. In this report, 
however, the lowest order gage functions are expressed in terms of their respective nat- 
ural parameters for purposes of simplicity. It will be shown that the gage functions of 
the thickness terms are of the same order of magnitude as those of some of the second- 
order lift terms. 

The boundary conditions (eqs. (6)) can be expressed in terms of the component 
potentials of the inner expansion (eq. (16)) to second order in 6 and sin o' as 
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(17) 


lllllllllll 


llllllllllllllllll 


9y 6 (x,r b> g) _ x dF (x) 

9r 6 dx 

d<p 6 (x,r^tO) _ d<p 6 (x,rrkn) _ 
90 90 


~\ 


±C w r 


9h(x,y) 


9x 


y 


9^l(x,r b ,0) 


9r 


= -sin 0 


90 1 (x,r,O) _ dcp^xf,*) 

90 90 


= -r 


(18) 


9g(x,y) 

ax 


j 


9^ 2 (x,r b ,0) 


ar 


= 0 


a0 2 (x,r,±O) a0 2 (x,r,±7r) 


dd 


d6 




j (x,r ? ± 0 ) 0 g(x ? y) j 90j(x,y^tO) ag(x,y) 

+ n “ 


ax 


ax 


ay 


ay 


(19) 




The potentials ^2 1’ anc * ^2 2 sa ti s fy homogeneous Neumann boundary 

conditions. 

The governing equations for the component potentials in the inner expansion (eq. (16)) 
are obtained after substitution of that expansion into equation (7) and collection of terms 
of the same magnitudes. It is found that the first-order potential <p^ and all the second- 
order potentials except <j ) 2 satisfy the two-dimensional Laplace equation in the cross- 
flow plane. For example, 


v 2^1 = 0 


( 20 ) 


where 


v2 = -2— + I -9- + -L -9=- 

2 9 f 2 r 9r f 2 9 0 2 


(21) 


The potential <^2 satisfies the Poisson equation 


*2*2 - ^ + + sin 0 + Fl^T cos e + + 


( 22 ) 
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The governing equations for the higher order potentials are discussed subsequently. It 
is interesting to note that the governing equations of the inner potentials are all parabolic 
since the determinant of the second- order coefficient matrix vanishes for each of these 
equations. (See ref. 14.) 

First-order potential .- The governing equation and boundary conditions for the first- 
order potential 0^ are given by equations (20) and (18), respectively. Solutions for this 
potential, which are traditionally termed slender-wing solutions, have been obtained pre- 
viously for the three wing-body configurations depicted in figure 5. 



(a) Attached leading-edge (b) Separated leading-edge (c) Attached leading-edge 

flow. flow. flow; swept trailing edge. 

Figure 5.- Wing-body configurations with known slender-wing cross- flow solutions. 

In this report flows in which the lower surface-wetting streamlines which approach 
the leading edge separate from the surface at the leading edge are termed separated 
leading-edge flows. Flows in which the surface-wetting streamlines remain attached to 
the surface are termed attached leading- edge flows. The solution for attached leading- 
edge flow past wing-body configurations with flat wings with straight trailing edges 
(fig. 5(a)) was obtained by Spreiter (ref. 15) and Ward (ref. 16). It should be noted that 
this solution is an extension of the basic slender-wing solution obtained by Jones (ref. 17). 
The solution for separated leading- edge flow past configurations of the same type (fig. 5(b)) 
was obtained by Wei, Levinsky, and Su (ref. 18) and is an extension of the work done by 
Mangier and Smith (ref. 19) and Smith (ref. 20) on separated leading- edge flow past slender 
flat-plate delta wings. The solution for attached leading-edge flow past wing-body config- 
urations with flat -plate wings with swept leading and trailing edges (fig. 5(c) was developed 
by Mirels (ref. 21) and Mangier (ref. 22). These solutions are reviewed in appendix A. 
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No slender- wing solution for separated leading-edge flow past a wing with a swept trailing 
edge is known to the author. 

As mentioned previously, the inner and outer expansions are to be matched in the 
region where the inner radial variable r is large. In appendix A it is shown that in the 
region beyond the wing where r » y 2 (x), the potential can be written as 

<Mx,r,0) = f(x)£%i (23) 

1 r 

For the configurations in figures 5(a) and 5(b), which have swept leading edges only, 


f(x) = \ y|(x) 


(24) 


where y 2 (x) is the function for the leading edge of the wing. It is shown in appendix A 
that this equation is valid to lowest order for both attached and separated leading- edge 
flow. For the configuration in figure 5(c), which has both swept leading and swept trailing 
edges, the function f satisfies the equation 


i\x) = S(x) 



y 2 (x) y 2 (x) 


(25) 


where K and E are complete elliptic integrals of the first and second kinds, respec- 
tively, with the argument 


k(x) = 




(26) 


The quantity S(x) is the same as the function S used by Mirels (see section 4 of 
ref. 21) and H used by Mangier (see section 3 of ref. 22) for wings with swept trailing 
edges. It is shown in appendix B that for configurations with twisted and cambered wings, 
the function f can be approximated as 


f(x) = y 2 (x) - | g'(x,z )|/y 2 (x) - V 2 dr] 


(27) 


if only the leading edge of the wing is swept. It is also shown that if both the leading and 
trailing edges of the wing are swept, the derivative f'(x) is approximately 


f'(x) = S(x) 



E(k) 

K(k) 


y 2 (x) y 2 (x) - i B q (x) 


(28) 
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where the modulus of K and E is given by equation (26). The function B 0 (x) is the 
same as the function B used by Klunker and Harder (ref. 23). The determination of this 
function is discussed in section 2 of reference 23. 


In the region very near the wing, the potential <p j can be expressed in the form 


0l(x,y,z) = m ± (x,y) - 


1 


ag(x,y) 

ax 


z 


(29) 


in order to satisfy boundary conditions (eqs. (18)). This form of the potential is used in 
the determination of the surface flow properties and some of the higher order terms. The 
upper and lower signs in equation (29) apply on the leeward and windward sides of the wing, 
respectively. 

Thickness potentials .- The terms in the inner expansion (eq. (16)) involving the 
thickness potentials cp 8 and <p 8 ^ are of second order in that the gage functions of 
these terms are 6^ and 6 2 log e (l/\v), respectively. Both potentials are governed by 
the two-dimensional Laplace equation in the cross-flow plane. The potential cp 8 sat- 
isfies boundary conditions (eqs. (17)) and the potential cp 8 j satisfies homogeneous 
Neumann boundary conditions. 

It has been shown by Stocker (ref. 24) that the solution for cp 8 is 


(p 8 (x,r,e) = g 6 (x) + ip(x,r,9) 


where 


m,r,9) = F(x) F'(x) 


2£w ri=y 2 ® 

\=r b (x) 


h'(x,7]) dr] 


+ ^ J _ ~ 2 h'(x,?7hlog e |/(r 2 + 7? 2 ) - 4r 2 r? 2 cos 2 9 

7 l =r b' x ' • 


+ log, 



4(x)\ 4r 2 r£(x) __ 2 


cos“ 9 / dr] 


(30) 


V* 


For r » y g (x) , equation (30) can be approximated as 


if/ = F e (x) f!(x) logp r 
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where F e (x), the nondimensional thickness distribution of the axisymmetric body with the 
same cross-sectional area, is written as 


F e (x) 


^F 2 (x) + 


4C W rv=y 2 m 
\=r b (x) 


h'(x,? 7 ) d7? 


(31) 


Since the governing equation and boundary conditions for cp^ ^ are homogeneous 
and are expressed in terms of the independent variables in the cross-flow plane, this 
potential can only depend on the independent variable x. Oswatitsch and Keune (ref. 1) 
have shown that if the leading term in the outer expansion depends on the thickness, as is 
the case for flows which are thickness dominated and flows for which the effects of lift 
and thickness are comparable, the potential (p g ^ is determined when the inner and outer 
expansions are matched. With equations (9) the thickness part of the inner expansion 
(eq. (16)) can be written in terms of outer variables for large values of the inner radial 
variable r as 


6 log e(i?r 6,1 (X) + 6 ‘ 


gg(x) + Fg(x) Fg(x) logg r 


= 62 l0 Se Gfe. 


0 6>1 (x) + F e (x) Fg(x) 


+ 5" 


g*(x) + F e (x) Fg(3f) logg T 


(32) 


It has already been shown that if the leading term in the outer expansion depends on thick- 
ness, the gage function of this term can be written as 

e x = v 2 = 6 2 

It follows that the leading term on the right-hand side of equation (32) must be zero so that 
< P 6}1 = -F e F; (33) 


Second- order lift potentials .- The potentials 02,1 an< ^ ^2,2 are governed by 
the two-dimensional Laplace equation in the cross-flow plane and satisfy homogeneous 
Neumann boundary conditions. Consequently, these potentials are functions of the inde 
pendent variable x only and can be written as 


^2,1 = S2,l® 

02,2 = S2 ; 2 


(34) 
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The potential 0 g is governed by equation (22), a Poisson equation, and satisfies boundary 
conditions (eqs. 19)). The particular solution for 0 2 was first presented in reference 8. 
This solution, which was derived for large values of the inner radial variable in terms of 
wind-oriented coordinates, is written in terms of the body-oriented coordinates used in 
this report as 


0 


p _ y + 1 

2 ~ ~ T ~ 


f'(x) f"(x) log| r - 



y + 1 

~T~ 


f 


"(x) 


cos 26 + 


m f'(x) 
2\ 2 ? 2 


(35) 


In reference 9 Cheng and Hafez present the solution for 02 * n t erms of wind- 
oriented coordinates in a form which is not restricted to large values of the inner radial 
variable. This solution is derived in detail in reference 25. The solution in terms of 
body-oriented Coordinates can be obtained in a similar manner. Let the potential 02 
be written as 


0 2 = g2 + 4>2a + ^2b + ^2c + ^2d 


(36) 


where the potentials 0 2a , 02b» and 02c are governed by the Poisson equations 


v|0 2a = (y + 1)^101 


(37) 


v 2 rf> - 2 Jmtk 4 . JL 

2^2b x 2y8r 9r d6 dd ) 


(38) 


V 2 0 2 ^ = sin 0 + 4 -^r cos 6 


,2\ dr 


r 86 


(39) 


and satisfy homogeneous Neumann boundary conditions. The potential 0 2c j is governed 
by the two-dimensional Laplace equation in the cross- flow plane and satisfies the bound- 
ary conditions 

^ 2 d(*’ ? b’ 0 ) _ n 
8? 


902d(x,r^°) 

86 


(40) 


90 2d (x,r,±7r) jr 

8<Pi(x,y,+0) ag(x,y) 1 9^l(x,y,+0) ag(x,y) 

86 " T x 2 

— 1 

CM 

h 

$8 


Consider the particular solution of 02 a and 0 2 iy 
the cross-flow plane in terms of inner variables be 


Let the complex variables in 
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— voi 0 


X = re 


and let the complex conjugate of X be designated as X*. Equations (37) and (38) can 
be written as 


and 


^2a _ 

ax ax’* 


= 4=fr; 

ax v 1 


(41) 


d_ 02b _ _ 1 _ 9< ^1 9< ftl 

ax ax* ~ a . 2 ax* 


(42) 


respectively. Cheng and Hafez (refs. 9 and 25) show that the general form of the partic- 
ular solutions to equations (41) and (42) are 




= r f i Ai’ dx i' dx *[^ (x ’ x ' x * ) ]' 


(43) 


^2b ~ ~2 ^1^1 


(44) 


More detailed forms of these solutions are given in references 9 and 25 and in appendix C 
of this report. 

It should be noted that the potential 02 c a PPears in the present treatment and not 
in that of references 9 and 25 because body-oriented coordinates are used in the former 
case and wind-oriented coordinates in the latter. Equation (39) can be written in terms 
of Cartesian coordinates as 


a2 ^2c _ 2 8 ^1 

" 2 az 2 > 2 3z 


dy“ V 

The particular solution to equation (45) is 


(45) 


A 


(46) 


,P 

^2a 


It has been shown by Cheng and Hafez (ref. 9) that the particular solutions 

P 

and 02b do n °t satisfy the appropriate boundary conditions. It is shown in appendix C 
that the particular solution 02 c does not satisfy the boundary conditions either. In 
order to enforce the boundary conditions, it is necessary to include the complementary 
In reference 25 Cheng and Hafez derive the complementary solutions 02 a 


solutions. 
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and 02b subject to the assumption that the presence of the body can be ignored. This 
assumption can be made since the body radius is, in general, small compared with the 
wing span so that the terms depending on the body radius are of higher order. The gen- 
eral form of the potential cf > 2a is written as 


log e 


|(s - X)(s - X*) ds 


(47) 



p C 

Similar expressions pertain for the potentials $ 2 b anc * ^2c' Detailed expressions for 
these solutions are given in appendix C. It should be noted that it is shown in reference 25 
that in order for the solution 0 2a to exist, the derivative <p ^ must be differentiable 
in y at all points on the wing. This condition is not met, in general, by attached-flow 
solutions for (j>y The solution for 0 2a derived i n references 9 and 25 was restricted 
to attached flow past wings with leading edges which are drooped in such a manner that 
the differentiability condition is met. Consequently, there is at most one angle of attack 
for a given wing at which a solution can be obtained. In appendix C it is shown that the 
solution presented in references 9 and 25 is also valid for separated leading- edge flow at 
arbitrary angles of attack. 

The potential 0 2c j> which is governed by the two-dimensional Laplace equation and 
satisfies boundary conditions (eqs. (40)), is written as 


^2d = " 


4n\‘ 


z 


=y 2 (x) 

=-y 2 (x) 


m'(x,s) g* (x,s) + 


1 8m(x,s) 9g(x,s) 


ay 


ay 


log e j(s - X)(s - X*)j ds (48) 


where 

m(x,y) = ijm + (x,y) - m_(x,y)j (49) 

The presence of the body has been ignored as before. It should be noted that the poten- 
tial 0 2( j does not appear explicitly in the treatments of references 9 and 25 because the 
lift potential is not resolved into first- and second-order parts in those references. 

From the results of appendix C, it is seen that for r » y 2 (x)» the potential </> 2 is 
the sum of the potentials <p ^ and <p^ where is given by equation (35) and 

is written as 
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aC „ f(x) f (x) cos 26 , „ ~ 

0 2 “ g 2 (x) - — ^2 ~2 + H(x) l0g e r 


4A 4 


( 50 ) 


where the coefficient H(x) is 


H (x) = -Jr J y'(x,s) P.V. j* ^ + log 1 1 - s|)dt ds 

1 b7r 9x J s=-y„x J t=-v 0 x V 1 s e v 


-y 2 (x) 


t=y 2 (x) t/ 
y 2 (x) 


1 ps=y„(x) 

-i-s \ m(x,s)ds — ^ » 

! 7 rA ZJ fi=-v ^ 27T\^ J s=-y 0 (x) 


27 T\* w s=-y 2 (x) 


i_ f s= y 2 (x) 

2 ( 


m 


'(*,«) e'M.-!,- !- s) 8g( ^’ s> 


9y ay 


ds 


1 ns=y 0 (x) ( r -T) 

\ \m(x,s) g"(x,s) - m'(x,s)ll - g'(x,s) > ds 

2ir \ 2 J s=-y 2 (x) < J 


(51) 


The function y(x,y) is the wing vorticity at the point x,y. 

The relationship between the radial stretching parameter v and the angle of 
attack a for flows for which the effects of lift are either dominant or comparable with 
the effects of thickness can be determined from a study of the second-order lift poten- 
tials. It can be shown from equations (9), (34), (35), and (50) that the sum of the second- 
order lift terms in the inner expansion can be written in terms of outer variables for 
large values of the inner radial variable r as 


log efe)* 2,2 ( *> + log e(r^2,l 


(x) + 0 2 (x,r,0)| 


= A 4 sin^ a log| , J_ 

e 


g 2,2 


(x) + 


y + 1 


f’(x) f"<x) 


n2olog efe 


X 


g 2,l 


(x) + H(K) + - * 1 f’ft) f"(x) log e r 


+ A. 4 sin 2 o |g 2 (x) + H(x) log e r 


- ^ - f’(x) f"(x) log 2 r - | f'(x)j 


y + 1 


f"(x) 


cos 20) + 


(52) 


If the effects of lift are to influence the first term in the outer expansion, the gage func- 
tion of the leading nonvanishing term on the right-hand side of equation (52) must be equal 
to ej. It has been shown that the parameter v can be related to by equation (13). 

In reference 8 it was assumed that none of the terms on the right-hand side of equation (52) 
vanish. This assumption led to the conclusion that a was related to \ and v by the 
equation 
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sin o' = 


at/ 

l0 Se h 


where a is an order- one constant. In addition, it necessarily followed that if the effects 
of lift and thickness were comparable, 5 was related to X and v by the equation 


6 2 = 




log, 


e \u 


where /i is an order-one constant, rather than by equation (14). In reference 9 Cheng 
and Hafez present a second solution in which it is assumed that the first term in equa- 
tion (52) vanishes identically. Consequently, the function g 2 2 ( x ) satisfies the equation 


^ 2,2 ~ 


y+J. f . f „ 


(53) 


and a can be related to \ and v by the equation 


sin oi ~ 


ay 


V log e h 


(54) 


It also follows that v and 6 can be related by equation (14) for the case where the 
effects of lift and thickness are comparable in the outer region. It is felt that this second 
solution is physically more realistic and, as a result, is the one used in this report. 

Higher order potentials .- In the strictest sense, the acceptance of equations (53) 
and (54) should be contingent on the proof that these equations, when transformed to outer 
variables with equations (9), do not cause the higher order terms in the inner expansion 
for large values of the inner radial variable to appear to be of lower order than ej. It 
is shown in appendix D that the sum of the transformed third-order terms in the inner 
expansion is not, in fact, of lower order than e^. The same relationship can be shown 
to hold for the sums of the fourth and higher order terms of the inner expansion. Con- 
sequently, it can be concluded that equations (53) and (54) are acceptable and that the inner 
expansion for large values of r can be written as 

-J> = \ sin a - 6 2 log e (i;)F e (x) F^x) + 6 2 g 6 (x) + F e (x) F^x) log e r - \ 4 sin 2 a log 2 f'(x) f"(x) 


+ \ 4 sin 2 a log e ^) g 2;1 (*) + A- 4 sin 2 o/g 2 (x) + H(x) log e r + 1 ~ f'(x) f"(x) log 2 r 


4A/ 




1 _ y + 1 


f (X) 


cos 20) + 


(55) 
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Outer Expansion 

Equation (10) gives the form of the outer expansion. It has already been shown that 
the gage function of the leading term is related to the radial stretching parameter v (the 
reciprocal of the nondimensional radial length scale in the outer region) by equation (13). 
The higher order gage functions in the outer expansion can be determined from an obser- 
vation of the inner expansion written in terms of outer variables for large values of the 
inner radial variable r. From equation (55), it is seen that this expression is written as 


$ = i^ 2 <'g 6 (x) + g 2 > i(x) + a 2 H(x) + 


F e (x) Fg(x) + a 2 f'(x) f"(x) 


log e r I 


ija f(x) sin 0 ^ 2 a 2 


+ ' V ^, (g 2 fe) + H(x) log e r + ~ 4 — f'(x) f"(x) log 2 r 


log 7 “ f l0ge *- v 

e \v 


+ if’(x) 


L _ ILll f ”(x) 

2 4 1 w 


cos 20 ) + 


(56) 


It is seen from equation (56) that eg and eg are given by the equations 


e 2 = 


i/^a 


log efe 


(57) 


e 3 " 


A 2 

l °Zeh 


(58) 


The governing equation for the potential #4 determined by the substitution of 
equation (13) into equation (11) is 

K $4 + V 2 $4 = (y + l)^^ (59) 

When expansion (eq. (10)) is substituted into equation (7) and the gage functions ej, eg, 
and eg are specified by equations (13), (57), and (58), it is found that the potentials $g 
and #3 are governed by the equations 

k *2 + ^2*2 = (r + D (* 1 * 2 ) + sin 0 + | ^ cos 0 ) (60) 
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K *3 + V 2*3 - (V + 1) + *2*2_ + si " « + f-lT “ s ») * J2 *1 

As pointed out by Messiter (ref. 4), equations of this type can be solved iteratively for 
r « 1. For example, equation (59) is written in the form 

= -K$ J + (y + 

First, the pertinent complementary solution for r « 1 is obtained as 
$1 = G 1 (x) + A(x) log e r 


where the functions Gj(x) and A(x) are not yet known. A higher order solution is 
obtained after substitution of equation (63) in the right-hand side of equation (62) and inte- 
gration of the resulting Poisson equation. The solution can then be written as 

= G-^x) + A(x) log e r + - A'a”? 2 log 2 r - A" + - * — - A'a” 


r 2 log e 7- If (GlA’)’ +G l G i r + ’ • * 


This procedure can be continued to higher orders if necessary. The pertinent solutions 
to equations (60) and (61) are found, in similar fashion, to be 


<f> 2 = b(x) + \A’ - f B” + BMG} - \ A'l / r log e r si 


sin 9 


+ - * ^ (b'A*) r logp r sin 6 + . . . 


> 3 = G 3 (x) + Cj(x) log e r + r j 1 B t B" log| r - " ~ 4 ~ B ") cos 20 + • • • 


where the functions G 3 (x), B(x), and C(x) are, at this point, arbitrary. 


Matching of Expansions 

Most of the unknown functions of x in the inner and outer expansions can be deter- 
mined from a matching of the inner expansions for large values of the inner radial variable 


L 




i 

4 


(r » ^ 2 ) with the outer expansion for small values of the outer radial variable (r « 1). 

As pointed out previously, there are two ways to accomplish the matching, by inspection 
and by the term-by-term procedure which has been formalized by Kaplun and Lager strom 
(ref. 12). 

In this section these two expansions will be matched initially by inspection because 
of the simplicity of this method. It will then be evident that the regions of validity of the 
two expansions do not overlap. In order to satisfy the traditional requirement (which is 
sufficient but not necessary) that the matching of two expansions be done in a region of 
common validity, an intermediate expansion with a region of validity which overlaps those 
of the inner and outer expansions will be introduced. This intermediate expansion will 
be matched individually with the inner and outer expansions. It will then be seen that the 
expressions obtained for the inner and outer expansions with the two matching procedures 
are the same. Because the intermediate- expansion procedure meets the sufficient (and 
hence necessary) requirements, it can be concluded that the inspection procedure meets 
the necessary requirements for this problem. 

Matching by inspection .- An examination of the inner expansion for r » yg and the 
outer expansion for r « 1 shows that the undetermined functions of x in those expan- 
sions can be related as 

g 6 + g 2l + H = Gj 

A = F e r e + a2f ' f " 

B = f 


S 2 = G 3 
C = H 

Consequently, the inner expansion for r » y^ and the outer expansion for r « 1 can 
be written as 


<p(x,r,6) = U„fU Sin a f(5t) - 6 2 log e ^F g (j) Fg(x) + + F e (x) Fg(x) log e r - \ 4 sin 2 a logf t 1 f'(je) f’ 


m 


A 4 sin 2 a log e j H(x) + A 4 sin 2 (y/g 3 (x) + H(x) log e r + y + 1 f'(x) f”(x) log 2 r - -Op 


> 2 4 


1 £"«) 


m 

2A 2 r 2 


cos 20 + 


m. f, < g > 

4\ 2 r 2 


(67) 
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and 


cp(x,r,6) = Uootf WGjOi) + 


F e (x) Fg(x) + - ~ a 2 f'(x) f"(x) 


logg r. 


2 

+ , = a — (f(%) + . . 




log exP 


sin 0 
r 


*) + + h w i°g e r 

log e xP l 


+ 2^f-if’«)f"Wiog|r-^l 


^- 2 4 if "w 


cos 20 + 


( 68 ) 


respectively. The functions and Gg cannot be determined with this method. 

Methods for determining these functions will be discussed subsequently. 


If the effects of lift and thickness are comparable in the outer region, the param- 
eter v in equations (67) and (68) can be equated to the equivalent thickness 6 and the 
angle of attack a can be related to v, and hence 6, by equation (54). If the effects of 
lift are dominant in the outer region, the parameter v is related to a by equation (54) 
(the order-one constant a can be given a value of one) and the thickness terms in equa- 
tions (67) and (68) can be considered to be small. 


It should be noted that regardless of whether the flow in the outer region is governed 
by thickness or lift or both, the lowest order term in the outer region is a source term 
(a term which varies as log e r). It can be seen from equation (68) that even when lift 
is dominant in the outer region, the dipole term ^the term which varies as = sin 0^ is of 
secondary importance compared with the source term. From equation (68) it is seen also 
that the source strength is of the form 



(69) 


when the effects of lift and thickness are comparable in the outer region and of the form 



(70) 


when the effects of lift are dominant. (For this case the constant a is equated to one, 
as discussed previously.) The quantity is the apparent nondimensional rate of change 
in the x-direction of the cross-sectional area of the lifting configuration. It is well known, 
of course, that when the effects of thickness are dominant in the outer region, the source 
strength is of the form 
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¥ ®a = ( F I)' 


(71) 


! 

s 


It should be noted that the function f(x) is the x- distribution of lift for the configuration 
and that the terms in equations (69) and (70) which depend on f provide the thickness- 
like flow due to lift which was anticipated intuitively and depicted in figure 2. 

Matching with intermediate expa nsion.- Kaplun and Lagerstrom (ref. 12) have estab- 
lished a procedure for matching expansions term by term. It is assumed in this proce- 
dure that matching can be performed only in regions where both expansions are valid. An 
examination of expansions (eqs. (67) and (68)) shows that these expansions cannot have a 
common region of validity since the leading term of the former contains a dipole and no 
source and the leading term of the latter contains a source and no dipole. Consequently, 
the procedure of reference 12 cannot be used to match expansions (eqs. (67) and (68)) 
directly. 

Although the use of the Kaplun and Lagerstrom procedure is not necessary in order 
for a matching of expansions to be valid, it is sufficient in order to establish validity. 
There are two ways in which this procedure can be applied to the present problem. One 
approach is to modify the definition of either the inner or the outer radial variable or both 
so that the expansions in terms of these variables have a common region of validity. The 
second approach is to introduce an intermediate variable in such a manner that the region 
of validity of the intermediate expansion overlaps the regions of validity of the inner and 
outer expansions. Actually, the distinction may be somewhat academic since the inter- 
mediate variable may be simply a generalization of the inner or outer variable. However, 
from a procedural point of view, it is the second approach which is used in this report. 

The intermediate variables x and ? are chosen to be 
x = x = x 

f = XyrlogN(r-) = rlogN(^) 

where 

0 < N < °° 

Note that r coincides with r for N = 0 and that the ratio r/r increases monoton- 
ically with N but is always small. The expansion in terms of r should coincide with 
the outer expansion for N = 0 and should match term by term with the inner expansion 
for large values of N. The intermediate expansion is written as 
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(73) 


<p = UooiJInxi^*®) + m 2 x 2 (x,M) + • • •] 


From an inspection of the inner expansion (eq. (67)) written in terms of ? and the gov- 
erning equation (eq. (7)) written in terms of the intermediate expansion (eq. (73)) and ? , 
it can be shown that the first eight gage functions are written in terms of v, X, and N 
as 
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For N > the orders of magnitude of these gage functions increase monotonically with 
2 1 

the index. For N < these gage functions are ordered as 


M 2 > M3 > Mi > M7 > Ms > M4 > M5 > M6 


(75) 
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From equations (7), (72), (73), and (74), it can be shown that the intermediate potentials Xj 
to X 5 are governed by the two-dimensional Laplace equation in the cross-flow plane and 
that the potentials Xg> Xq, and Xg are governed by the Poisson equations 
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(78) 


It should be noted that the governing equations of the intermediate potentials are all inde- 
pendent of the parameter N. Consequently, the solutions are independent of N also. 
Higher order gage functions and the governing equations for higher order potentials can 
be obtained in a straightforward manner. 

When the Laplace equations governing the potentials x^ to Xg and equations (76), 
(77), and (78) for x 6 > Xy, and Xg> respectively, are integrated, and the resulting inter- 
mediate expansion with N ^ ^ is matched term by term with inner expansion in the man- 
ner outlined by Van Dyke (ref. 13), it is found that the solutions for the intermediate poten- 
tials are written as 
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Note that the validity of these solutions does not depend on the parameter N and hence 
the relative orders of magnitude of the intermediate gage functions. It can be shown that 
for N < the gage functions are ordered as inequalities (eq. (75)) indicate, and the inter- 
mediate expansion can be matched term by term with outer expansion (eq. (68)). (See 
eq. (65) for additional terms in the outer potential $2*) For N = 0 the intermediate 
and outer expansions coincide. 

It has been shown that there is an intermediate expansion (eq. (73)) with gage func- 
tions and potentials given by equations (74) and (79), respectively, which has a region of 
validity that overlaps those of the inner and outer expansions. It has also been shown 
that the results for the inner and outer expansions obtained by use of the intermediate 
expansion are the same as those obtained by inspection. Consequently, it is concluded 
that the inspection procedure is adequate for the present problem. 


Solution Near Configuration Surface 

The solution near the surface of the configuration can now be determined to within 
an arbitrary additive function of x. This solution is written in terms of inner variables 
to second order in sin a and 6 as 



where is one of the basic lift potentials discussed in appendixes A and B, f is the 

related dipole strength distribution, F e is the equivalent thickness distribution given by 
equation (31), ^ is the near-field thickness function given by equation (30), H is the 

secondary second-order outer-region source- strength distribution given by equation (51), 
02* a j and $2 are the second-order particular solutions given by equations (Cll), 

(C15), and (C20), respectively, $2a’ ^2b’ anc * ^2c are secon d-order complementary 
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solutions given by equations (C13), (C18), and (C23), respectively, and is the solu- 

tion given by equation (48). It should be noted that the presence of the body has been 

P P P C C C 

ignored in the derivation of the expressions for $ 2 a> <P 2b> P2c> $ 2 a> P2b> P2c> 
and 02d* ^ s h° u ld also be noted that the function 6 ^ G^(x) + X^ sin^ a Gg(x) cannot be 
determined with the present method. The determination of this function is discussed in 
the next section. 


Equation (80) for the perturbation velocity potential can be used to determine the 
pressure coefficient at the configuration surface. It is shown in reference 4 that the 
pressure coefficient in terms of body- oriented cylindrical polar coordinates is 
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Determination of Additive Function 

The solution for transonic flow over a lifting slender configuration has been deter- 
mined to within an arbitrary additive function of distance along the axis. It has been 
pointed out previously that this function cannot be determined with the method of matched 
asymptotic expansions. Consequently, the solution for this function must be obtained with 
another method. 

It should be noted that the method of matched asymptotic expansions can be used to 
determine the solution for potential flow past a slender body to within an arbitrary function 
of distance along the axis also. For subsonic and supersonic slender-body flows, which 
are governed by a linear equation, this function can be determined with the integral equa- 
tion method. (See ref. 26.) Spreiter and Alksne (ref. 27) and Aoyama and Wu (ref. 28) 
present approximate methods for determining this function for transonic slender-body 
flows, which are governed by a nonlinear equation. However, these approximate methods 
are not, in general, as accurate as required. At present, solutions of sufficient accuracy 
can be obtained only with numerical methods which treat the full two- variable problem 
rather than the one-variable problem for the additive function. One such method for treat- 
ing the transonic slender-body problem is that of Bailey (ref. 29), which is based on the 
successive line overrelaxation technique developed by Murman and Cole (ref. 30). 

The approximate methods of Spreiter and Alksne (ref. 27) and Aoyama and Wu 
(ref. 28) could be generalized and applied to the present problem of transonic flow past a 
lifting slender configuration. However, the more accurate approach of solution by numer- 
ical integration is preferable. Three-dimensional numerical calculations of transonic 
flow about wings and wing-body configurations have been made with methods based on 
that of Murman and Cole. (See ref. 31, for example.) However, because of their three- 
dimensional nature, these calculations require large amounts of computing time and 
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computer storage. Rapid two-variable approximate methods for treating the present 
problem have been developed by Barnwell (ref. 32) and Cheng and Hafez (ref. 25). The 
approximations used in these methods are based on the theory presented in this report 
and reference 9. It should be noted that the method of Barnwell can be used in the angle- 
of-attack range of order . _ = (6 is the equivalent thickness ratio), and that the 

V°Se | 

method of Cheng and Hafez is restricted to the much smaller angle- of- attack range of 
order 6^. 


NONSLENDER -C ONFIGURATION PROBLEM 
Wings are described traditionally as being nonslender if the reduced- span — length 


ratio 


i - m£ 


|b j SL is of order one. In this section nonlifting nonslender transonic wing 
theory is reviewed, the scaling relationships for transonic flow about nonslender lifting 
configurations are derived, and it is shown that second- order nonslender wing theory 
is in agreement with sweep theory although first-order theory is not. The agreement 
between second- order theory and sweep theory holds for both nonlifting and lifting flows. 


Nonslender Nonlifting Wing Theory 

Messiter (ref. 4) and Cole and Messiter (ref. 11) treat the problem of transonic flow 
past a nonlifting nonslender wing. This problem is characterized by the parameters 
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(81) 


It can be shown that for transonic flow past nonslender wings, there is no inner 
region where the governing equation simplifies. As in slender configuration theory, the 

radial length scale in the outer region is of the form H/v , and the quantities u, e. , and 

2 1 
1 - M» can be related by equation (13). (The slenderness approximation was not made 

in obtaining this equation.) It can be seen from equations (13) and (81) that the semi- 
span b, which is the radial length scale near the configuration, is of the form 



Consequently, the radial length scales in the inner and outer regions are of the same order 
and, as a result, the governing equation and scaling relationships for the outer region apply 
all the way to the configuration surface. 
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The relationship between the wing thickness-length ratio and the parameters v, 
ej, and 1 - can be determined from the surface boundary condition, which is of the 
form 


w c (x,y,z w ) = /tgX 
U ro \l) 


and the expression for the perturbation velocity component w c in terms of the perturba- 
tion velocity potential, which is 


w. 


9 w TT „ 9^1 


8z 
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It follows from the boundary condition and equations (13) and (82) that v and t w /fL can 
be related as 


v = 



(83) 


The gage function and the free-stream Mach number 
t -yj/lt by the equations 

2 / v 2 / 3 

1 - M* ... ft w 

1 K \i ) 


Moo are given in terms of 


(84) 


It should be noted that the scaling relationships for transonic nonslender-wing theory 
are the same as those for transonic two-dimensional theory. 


Nonslender Lifting- Wing Theory 

Consider the problem where the angle of attack a is small but is much larger than 
the wing thickness-length ratio t w / L The boundary condition at the surface of the wing 
is of the form 


w c (x,y,z w ) 

Uoo 


= 0(sin a) 


(85) 


From equations (13), (82), and (85), it can be seen that v is related to the angle of attack 
by the equation 

v = (sin a) 1 / 3 (86) 
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and the gage function 


and the free- stream Mach number are given by the equations 


e l 

£i = i^mI,oW /3 (87) 

If the wing thickness-length ratio t and the angle of attack a are of the same order 
of magnitude, these quantities can be related by the equation 

sin a = a —■ (88) 


where a is an order-one constant, and the quantities v, e^, and can be evaluated 

with either equations (83) and (84) or (86) and (87). 


It was shown in the previous section that the scaling for the outer region applies all 
the way to the wing surface for nonslender wings and that, as a result, the wing semispan- 
length ratio X varies as \/v. Consequently, for lift-dominated and thickness- dominated 
nonslender flows, X is scaled as 

X = j- (C = 0(1)) (89a) 

(sin a ) l/ 3 

and 


X = 



(C = 0(1)) (89b) 


respectively. When sin a and t w /j£ are related by equation (88), both of these equa- 
tions apply. 


It can be shown from equations (7), (9), (10), (83), (84), and (88) that when the effects 
of lift and thickness are comparable, the outer potential for nonslender flow is gov- 
erned by equation (59). 


Approximate Governing Equation for Nonslender -Configuration Problem 

It has been observed by Lomax, Bailey, and Ballhaus (ref. 31) that equation (59) is 
not adequate for treating flow over swept wings in that the conservation form of this equa- 
tion does not lead to the same inviscid shock-jump condition for the pressure for flow past 
an infinite swept wing as obtained from simple sweep theory. In reference 31 it is shown 
that the proper shock-jump condition for an infinite swept wing can be obtained if several 
terms are added to the governing equation. It can be shown that all of these additional 
terms are in the second- order approximation to the governing equation for cp. It is prob- 
able that other second- order terms should be included since the terms added in refer- 
ence 31 were chosen for a special case. 



The governing equation for the second- order potential $2 in the outer expansion 
(eq. (10)) for cp for a nonslender wing has been derived by Messiter (ref. 4). This equa- 
tion can be written as 
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The gage function was found to be 
e 2 = 1/4 


(91) 


The value of the quantity v in equation (91) is given by equations (83), (86), or both. Let 
the velocity potential be approximated by the first two terms in the outer expansion 
(eq. (10)) so that 


<P* e 1*1 + e 2*2 
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It can be shown from equations (13), (59), (84), (90), (91), and (92) that the governing equa- 
tion for the second- order approximation to <p is 
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It is clear that the terms on the left and the first term on the right are those obtained from 
the first-order equation. The second-order terms which were added in reference 31 are 
the terms on the right which involve derivatives with respect to y. 


COMPARISON OF SLENDER- AND NONSLENDER- WING THEORY 


In this section the Mach number range and flow- field structure for slender- wing 
and nonslender- wing flows are compared, and the Mach number range and flow-field struc- 
ture for flows which exhibit both slender-wing and nonslender-wing characteristics are 
determined. 
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Mach Number Range and Structure of Slender- Wing Flows 

For purposes of comparison, configurations composed of wings alone are considered. 
It has been shown that if slender-wing theory is applicable, the velocity potential in the 
outer region can be written to lowest order as 

<p = Uoo-W 2 $ 1 (x,y,^) (94) 

It can be shown from equations (1) and (3) that the equivalent- body thickness 6 can be 
written as 

62 = C^ X ^ (C w =0(1)) (95) 

where, for slender-wing theory, 

* = j = 0 ( 1 ) 

Since the parameter \ can be associated with the aspect ratio, the wings under consid- 
eration have aspect ratios of order one. It can be shown from equations (94) and (95) that 
the velocity perturbations in the outer region can be written as 



These perturbations are found in a region of the cross-flow plane bounded on the outside 
by a radius r of the order 



and on the inside by a radius which satisfies the inequality 
r » b = \i 
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Inside this region the perturbations are obtained from the inner expansion and, as a result, 
are larger. The governing equation in the bounded outer region where equations (96) apply 
is of mixed elliptic-hyperbolic type. In the inner region the governing equation is para- 
bolic and in the far field beyond the bounded outer region the equation is elliptic if 
is less than 1 and hyperbolic if is greater than 1. The Mach number range in which 

slender-wing theory applies is 


Moo =1 + KS 2 =1 + o(-2 


(98) 


Mach Number Range and Structure of Nonslender-Wing Flows 
It has been shown that the velocity potential for nonslender-wing theory can be writ- 
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ten as 
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and that X is of the form 
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Consequently, nonslender-wing theory applies to large-aspect- ratio wings. From equa- 
tion (99) it can be shown that the velocity perturbations are of the form 
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These perturbations are applicable in a region which extends above and beneath the wing 
to a distance of the order 


f-° 
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-1/3 


(102) 
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The governing equation is of mixed elliptic-hyperbolic type in this region and of elliptic 
or hyperbolic type in the far field. The Mach number range in which nonslender-wing 
theory applies is 


Moo = 1 + K 




2/3 


1 + O 



(103) 


It can be seen from equations (96) and (101) that for wings with a given thickness 
ratio, the velocity perturbations obtained from nonslender-wing theory, which applies to 
wings of large aspect ratio, are larger than those obtained from slender- wing theory, 
which applies to wings with aspect ratios of order one or less. From equations (97) 
and (102) it can be determined that the radial extent of the slender-wing perturbations is 
larger, and from equations (98) and (103) it is seen that the Mach number range in which 
the nonslender-wing perturbations occur is larger. 


Flows Exhibiting Both Slender-Wing and Nonslender-Wing Characteristics 

It can be shown that slender- and nonslender-wing theory are both applicable to flow 
fields about large aspect ratio swept wings if M^ is close enough to 1. Such a wing is 
shown in figure 6. The quantities L, b, and £ are the total length, semispan, and 
chord of the wing. It is assumed that 

A = £ = 0(1) 

x ’ it » 1 

If both theories apply to the same flow field, nonslender-wing theory applies close to the 
wing where the characteristic length scale in the x-direction is the chord £, and slender- 
wing theory applies at large distances from the configuration where the length scale in 
the x-direction is the total length L. 

In the region close to the wing where nonslender-wing theory applies, the velocity 
potential is given by equation (99) and the perturbation velocities satisfy equations (101). 
These equations pertain in the Mach number range given by equation (103) and in a region 
which extends above and beneath the wing a distance of the order of the length scale 
f(t w /f)" : *' // ^. This region of the cross-flow plane is depicted in the cross section in fig- 
ure 6. The governing equations in this region are of mixed elliptic- hyperbolic type. The 
aspect ratio of the wing is given by equation (100). It follows from equation (100) that the 
chord-length and thickness-length ratios are 
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From figure 6(a) it can be seen that the equivalent thickness-length ratio is given by the 
equation 
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Slender- wing theory applies if the Mach number is in the range 
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The velocity potential in the outer region (the region at large radial distances from the 
body) is written as 
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where, in this case, the outer variables x, y, and z are written as 


x = x 
X L 


3/2/0 5/6 


= L 


C U/ 

It follows that the velocity perturbations in the outer region are of the order 
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(a) Wing planform. 
Elliptic or hyperbolic far field 



Figure 6.- Large-aspect-ratio wing to which both slender- and nonslender-wing theory 

apply for = 1 + ofe) 5/3 ]. 



As shown in figure 6(b), the outer region is bounded on the outside by a radius of the order 



and on the inside by a radius of the order 
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According to slender- wing theory, the governing equation is of mixed and parabolic type 
in the outer and inner regions, respectively, and of either elliptic or hyperbolic type in 
the far field beyond the outer region. It can be seen from equation (102) and inequal- 
ity (105) that the mixed flow region near the wing predicted by nonslender-wing theory is 
embedded in the parabolic inner region predicted by slender-wing theory if the Mach num- 
ber falls within the range given by equation (104). 


CONCLUDING REMARKS 


A solution has been obtained for the perturbation velocity potential for transonic 
flow past lifting configurations with span-length ratios of order one. This solution per- 
tains to configurations for which the reduced span-length ratio (the product of the span- 

length ratio and the quantity |/|l - m£( where M^ is the free- stream Mach number) 
is small. The angles of attack which are considered are small but are large enough to 
insure that the effects of lift are either dominant or comparable to the effects of thickness 
in the outer region. The analysis was performed with the method of matched asymptotic 
expansions. 

It is shown that the lowest order effect of lift in the outer region is in the form of 
a source, and that the doublet effect is of secondary importance in this region. This find- 
ing pertains both when the effects of lift are comparable with those of thickness and when 
the effects of lift are dominant. As a result of the source flow due to lift, streamlines 
are deflected outward more than they would be by thickness effects alone. 

A short study is made of the flow about lifting nonslender configurations (configura- 
tions with reduced span- length ratios of order one). The order of magnitude of the veloc- 
ity perturbations is established, and it is shown that second- order nonslender-wing theory 
is in agreement with sweep theory. The results of slender- and nonslender-wing theory 
are compared, and the Mach number range and general flow-field structure of flows exhib- 
iting both slender and nonslender characteristics are determined. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 11, 1975. 
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APPENDIX A 


EXISTING CROSS- FLOW SOLUTIONS FOR WING- BODY CONFIGURATIONS 

As discussed in the text, solutions have been obtained for the cross-flow velocity 
potentials for the three wing-body configurations depicted in figure 5. The purpose of 
this appendix is to review these solutions and to determine the asymptotic forms at large 
distances from the configurations. 

Attached Leading-Edge Flow Past Configurations With Swept Leading Edges 

General solution .- A cross section of the wing-body configuration is shown in fig- 
ure 7. The solution for attached leading-edge flow past this configuration was obtained 

z 



t 


Figure 7.- Cross section of configuration composed of circular body 

and flat- plate wing. 

by Spreiter (ref. 15) and Ward (ref. 16). Let the complex variable in the cross-flow plane 
be 


X = y + iz 


The complex perturbation potential for this problem is written as 
W(x,X) = i|x - Y(x,X)| 

where 


Y(x,X) = 




(Al) 


(A2) 
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APPENDIX A - Continued 


The cross- flow perturbation velocity potential is 
0 1 (x,r,0) = Rejw(x,X)| 
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where the relation 

X = re i& 
has been used. 

Solution far from configuration .- For large values of r, equation (Al) can be writ- 
ten as 
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Thus, the velocity potential for large values of r can be written as 

r»“ 


= 2*2 1 + df 


b ' sin 0 


(A3) 


For practical purposes the quantity [? b/ /y 2 ) can be neglected since ? b and y 2 are 
generally of orders 6 and 1, respectively. Consequently, the velocity potential for 
large values of r is 



APPENDIX A — Continued 

Separated Leading-Edge Flow Past Configurations With Swept Leading Edges 

General solution .- Consider configurations composed of circular bodies and flat- 
plate wings. The cross- flow solution for nonconi cal configurations of this type was 
obtained by Wei, Levinsky, and Su (ref. 18) by use of the method developed by Mangier 
and Smith (ref. 19) and Smith (ref. 20) for conical flat-plate wings. The physical cross- 
flow plane and the transform plane are shown in figure 8. The outer part of the vortex 



t 

(a) Physical plane. 

Y -plane 



t 


(b) Transform plane. 

Figure 8.- Model for separated leading-edge flow. 
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APPENDIX A - Continued 


is approximated with a segmented vortex sheet which emerges from the wing leading edge. 
The inner part of the vortex sheet is modeled with a vortex core and a feeding sheet which 
connects the core to the vortex sheet. 

The complex perturbation potential for this model is written as 
W(x,X) = i(x - Y(x,X) - ± r v (x)log e 


Y(x,X) - y(x,X v ) 
Y (x,X) + Y (x,Xy j 


. o= oi 


r — max 

sL 




Y(x,X) - Y|,X(q]| 
e \Y(x,X) + Y[X,X*(a)] 


da 


(A5) 


where the asterisk denotes the complex conjugate, X v and X(a) are the positions of 
the vortex core and the centers of the vortex- sheet elements, T y is the strength of the 
vortex core, and A 0 is the potential jump across the vortex sheet. The values of these 
quantities are determined from the leading-edge Kutta condition, the condition that the 
pressure and normal velocity are continuous across the vortex-sheet elements, and the 
assumption that no force is exerted on the system composed of the vortex core and feeding 
sheet. 

Solution far from configuration .- It can be shown that for angles of attack such that 
sin a « 1 the term proportional to T v and the integral term in equation (A5) are of 
order (sina) c where c £ 1 so that equations (A3) and (A4) apply to separated as well 
as unseparated leading-edge flow. Let Y(X,X V ) and Y(x,Xy) be written as Y v and 
Yy, respectively. In the region where r and hence Y(x,X) are large, the logarithmic 
expression for the vortex core in equation (A5) can be written as 


l°g e 


' Y - Yt 


V Y + Y* 



- log e 




Assume that the vortex location is 


(A6) 


Xy = y 2 - d v + iz v 

where d v and z v are shown in figure 8. It is reasonable to assume that for small 
angles of attack, the distance between the vortex core and the wing is of the form 

z v = 0(sin a) 
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In fact, Brown and Michael (ref. 33) show that this is the case for their simplified model 
for conical flow past a delta wing. For d v , there are two cases of interest. If the vor - 
tex core lies near the leading edge of the wing, then d v is of order (sin a) c where 
c g 1 (for the Brown and Michael model for conical flow, c = If the vortex core lies 
inboard of the tip, d v is of order one. 

First consider the case where 

d v = o(sin c a) (A7) 


After some algebraic manipulation, it can be shown that Y v is 


where 


Y v = 


r4 


v. i - !} Va* ♦ 2? a** 


/3 = ir + tan-lf^^' 


0 = tan-1 52 


(3y>0) 

(3v<o) 


(A8) 


(A 9a) 
(A 9b) 


and that 


0 

cosf = 


^/dy + Zy “ dy 


# 


+ zt 


(A10) 


It follows from equations (A8) and (A10) that for large values of r, equation (A6) can be 
written as 


r y - y v , 
1o My7y*J " 2 




, y 2 | 1 -^rv + ^- d v)V- = v-o 
Hi 


e“ i6f e~ id 


l-$ 


(sin a) 


In a similar fashion the terms log e 


Y - Y(ct) 


can be shown to be of this order of mag- 


[Y - Y*(a)| 

nitude also. It can be seen from reference 18 that the leading- edge Kutta condition is 
written as 
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Myv). r%-d W !fM d „ = , 

{ Y v | 2 *^=0 daL J |y(ct)| 2 


(All) 


From equations (A7), (A8), and (A10), it can be determined that 
Re(Y v ) 


= O 


(sin a) 


i 3c 
2 


(A12) 


Re[Y(o)] 

A similar expression applies for — r-. Since the vortex- core term and the integral 


|Y(o)|‘ 


term in equation (All) are generally of the same order of magnitude and do not cancel, it 
follows that 


r v = o 


(sin a) 


3 c i 
2 


1^=0 

da 


(sin a) 


3c i 
2 


This order- of- magnitude estimate for r v is consistent with the delta- wing results of 
Brown and Michael (ref. 33) that c = ^ and r y is of order 1. Consequently, the 
vortex-core term and the integral term in equation (A5) are of the order (sina) c and 
can be neglected in the first approximation. 

Consider the case in which 


d v = 0(1) 


and let 


h 


- d„ 


After considerable manipulation, it is found that 


= l/yf-fV 


i- 


b_ e i/3/2 


y2ar2 

y 2 y v 


(A13) 


where 


/3 = ir + tan“l 




" 2z v y 2\ y v ‘b 


y v ( y i - y v)( y l y v - *b) 


(A14) 
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and that 


p 

C0S 2 = 


W - ? b 


H y l ' wPv ' ?4 


(A15) 


From equations (A6), (A13), and (A15), it can be shown that in 


log. 


2 (^ - 


a -i0 a -id 


\Y + Y, 




the first approximation, 


0(sin o) 


Similarly, the quantities log 


Y - Y(or) 


Y + Y*(cr) 


can also be shown to be of this order. The 


magnitude of the core strength T y and the quantities -dA0(o)/dcr cannot be estimated 
from equation (All) in the manner used previously because the orders of magnitude of 
the distances from the vortex core and vortex- sheet segments to the wing tip vary and 
the orders of magnitude of the coefficients of the quantities -dA<£/dcr and r v in equa- 
tion (All) depend on these distances. However, it is reasonable to assume that the core 
strength T v and the vortex-sheet segment strengths -dA0/dcr are of order 1 or 
smaller. Consequently, the vortex- core term in equation (A5) is of order sin a or 
smaller, and the contributions of the various vortex-sheet segments to the integral term 
vary from order sin a or smaller near the core to order (sin ct) c near the tip. These 
contributions can be neglected in the first approximation. 


It is concluded that for separated leading- edge flow, the velocity potential for large 
values of r is given to lowest order by equation (A3) or (A4), which were derived for 
attached leading-edge flow. It has been shown that this conclusion applies both when the 
vortex core is located near the wing tip and when it is located well inboard. 


Attached Leading- Edge Flow Past Configurations With Swept 

Leading and Trailing Edges 

General solution .- The derivative 
Mirels (ref. 21) and Mangier (reL 22). 
obtained with the transformation 

r^ 

z = X - a? (A16) 


for this problem has been derived by 
The physical plane and the transform plane 
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are shown in figure 9. 


In the transform plane the derivative 


9 2 W 

axaz 


is written as 


9 2 W 

9x9Z 


-i?7 2 ?72 S(x) 




(A17) 


where 


*1 = y i - r 


Vo = yo-^r 



t 

(a) Physical plane. 


h Z plane 

i 2if b 


o 





o- 

^1 



v 


A -2i? b 

t 


(b) Transform plane. 

Figure 9.- Cross-flow planes for attached leading-edge flow past configuration with 

swept trailing edge. 
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The function S(x) is the function Mirels (ref. 21) developed for enforcing the Kutta con- 
dition at the swept trailing edge, and K and E are complete elliptic integrals of the 
first and second kinds with the modulus 



It should be noted that Mirels’ S- function is identical to the function H used by Mangier 
(ref. 22). 

Solution far from configuration .- The complex velocity in the physical plane is 
related to that in the transform plane by the equation 


3W _ aw BZ 

3X az 9X 



It follows that the equation for 


a 2 W 

ax ax 


is 


a 2 w 
ax ax 



2 aw ~ 


x 


2 az 


r b r b 


(A19) 


At large distances from the wing and body in the transform plane, equation (A17) can be 
written as 


£& = -iT7 0 77;4rS(x)fl - § 


axaz 


'2 '2 z 2 


K 


and hence the derivative aW/aZ is approximately 


aw = , 
az 




= 1 


11 + V 


Thus equation (A19) can be written as 


dfW = _ _i_ 


ax9X 


X* 


S(X) l 1 ■ ijV2 ■ 2F b F b 
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so that the derivative aW/ax is written as 



Therefore, for large values of r the derivative a^j/ax can be written as 



(A 20) 


where the modulus of E and K is given by equation (A18). The quantities r^ and y 2 
are of orders 6 and 1, respectively, and the quantity yj varies between orders 5 
and 1. Consequently, equation (A20) can be written to lowest order as 




where the modulus of E and K is 



(A22) 




APPENDIX B 


FAR- FIELD APPROXIMATION TO CROSS- FLOW SOLUTIONS FOR 
CONFIGURATIONS WITH TWISTED AND CAMBERED WINGS 


The purpose of this appendix is to develop the far-field approximation to the cross- 
flow solution for configurations with twisted and cambered wings. The wings which are 
treated include those with only swept leading edges and those with both swept leading and 
swept trailing edges. As in the case of flat wings, the approximation for the cross-flow 
potential (p j can be determined analytically for wings with swept leading edges only, but 
only the derivative 8$jyax can be determined analytically for wings with both swept 
leading and swept trailing edges. The velocity potential is governed by equation (20), the 
two-dimensional Laplace equation in the cross-flow plane, and satisfies the boundary con- 
ditions (eqs. (18)). Let the potential be written as the sum 


^1 ~ ^la + ^lb 

where the potentials and $j b satisfy the boundary conditions 


8 ? 


= -sin 0 


(Bl) 


(B2a) 


and 


801a(x,r,°) _ 80 la (x,r,jr) 
80 80 


(B2b) 


80ib(g,r b ,g) 

8r 


(B3a) 


8</> lb (x,r,0) 

80 


80 lb (x,r,7r) 

80 


r g’(x,y) 


(B3b) 


respectively. Equations (B2b) and (B3b) apply on the interval y 2 = r ^ ? b if the wing has 
no trailing edge and on the interval y 2 § r g yj if the wing has both leading and trailing 
edges. If the wing has a trailing edge, equations (B2) and (B3) are supplemented by the 
equations 


80 la (x,r,O) 

8K 


80 la (x,?,Jr) 

8x 


(B4) 
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and 

90 lb (x,r,O) = d<t> lb &,r,n) = q 
9x ax 

respectively, on the interval yj > r i r^. 

Configurations With Swept Leading Edges 

The solution for the potential <j !>^ a which satisfies boundary conditions (eqs. (B2)) 
for attached leading-edge flow past configurations with unswept trailing edges is derived 
in the first part of appendix A. The far-field approximation for this potential, which is 
given to lowest order by equation (A4), is also derived. It should be noted that since 
equation (A4) does not contain r b , the use of the lowest order form is equivalent to 
neglecting the presence of the body. In a similar fashion the far-field approximation for 
the potential which satisfies boundary conditions (eqs. (B3)) can be approximated 

with the solution of Klunker and Harder (ref. 23) for attached leading-edge flow past 
twisted and cambered wings. Let the complex variable in the cross-flow plane for the 
wing alone be X, and let the boundary condition for the wing be 


w c (x,y,0) = g'(x,y) 
-y 2 (x)sy sy 2 (x) 


From equation 1.3 of reference 23 it is seen that the complex velocity in the cross-flow 
plane of the wing is 


9W(x,X) 


ps=y 2 (x) g’(x,s)|/y2(5c) - s 2 
\ — — ds 


_ 


8X 


771 /x 2 -y|(x) J s=-y 2 (S0 


X - s 


(B6) 


At large distances from the configuration equation (B6) can be written as 


MgS . _2i_ f s=y ~2 R) g . (x ,s)&7^ ds 
9X nX ZJ s=0 v 1 


(B7) 


The velocity potential (p^ for large values of r can be obtained from equation (B7), 
and the potential can be obtained from equation (A4). Consequently, at large dis- 
tances from the configuration, the potential is written as 
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$ 1 (x,r,e) = 


2 y 2' 'll? 2 g ' ( *’ s) l/ i 1" = 2 dS 


sin 0 


(B8) 


Configurations With Swept Leading and Trailing Edges 

The solution for the derivative 9 ^la/ 9 ^ °f the potential which satisfies boundary 
conditions (eqs. (B2) and (B4)) for attached leading-edge flow past configurations with 
swept leading and trailing edges is derived in the last part of appendix A. The solution 
for the derivative of the potential which satisfies boundary conditions (eqs. (B3) 

and (B5)) can be approximated with the solution of Klunker and Harder (ref. 23) for flow 
past twisted and cambered wings. The boundary conditions for the wing treated in refer- 
ence 23 are written in terms of the notation of this report as 

w c (x,?,0) = g’(x,y) 



u(x,?,0) = 0 


-yj ^ y g y 1 


Let the complex variable in the wing cross-flow plane be X. From equation 1.8 of ref- 
erence 23, it is seen that the x derivative of the complex velocity in the cross-flow plane 
of the wing is 


9 2 W(x,X) 

axax 



’(x,s)^y 2 - s 2 


3/2 




2 . ?2 


X - s 


ds 


\3/2 


1 


r: 


s=y 2 g"(x,s)(V| - s 2 ) |/s 2 - yf 


s=?i 


j_ 

X - s 


ds + A 0 (x) + B 0 (S) X“ 


(B9) 


where A 0 (x) and B Q (x) are the functions A and B used in reference 23. 
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At large radial distances from the configuration, equation (B9) can be written as 


9 2 W . i B off) 
9x9X ^ 


(BIO) 


The x derivative of 0-^ can be obtained from equation (BIO), and the derivative of 
0j a can be obtained to lowest order from equation (A20). Consequently, the derivative 
of 0^ for large values of F can be written as 


801 

ax 



sin 6 
r 


where the modulus of the complete elliptic integrals E and K is given by equa- 
tion (A22). A method for determining the function B 0 (x) is discussed in section 2 of 
reference 23. 
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SOLUTIONS FOR SECOND- ORDER LIFT POTENTIALS 0 2a > 02b’ AND 02c 

Particular solutions for the potentials 0 2 a > 02b’ and 02c are derived in very 
general form and are given by equations (43), (44), and (46), respectively. The purpose 
of this appendix is to evaluate these general forms at the configuration surface and to 
determine the complementary solutions necessary to enforce the boundary conditions. 

Cheng and Hafez (ref. 25) have shown that in order for the solution for 0 2a 
exist, the potential 4 > j must satisfy a Kutta condition at the leading edge of the wing. In 
reference 25 solutions are obtained for the potentials 0 2a and $ 21o for configurations 
with wings which are cambered so that this condition is satisfied. In this appendix it is 
shown that the solutions derived in reference 25 also apply approximately to flow fields 
where the leading- edge Kutta condition is met by means of leading- edge separation. 

As stated in the text, the boundary conditions for the potentials 02a » 02b’ anc * 02c 
are all homogeneous. In this appendix the body- surface boundary condition is ignored as 
it was in reference 25. This approximation is permissible since the body surface r^ is 
of order 6. 


Cauchy Representation of <p j and for the Leading-Edge 

Separation Model 

The Cauchy representation is used in this appendix as it was in reference 25. The 
leading-edge separation model which is employed is that used by Mangier and Smith 
(ref. 19), Smith (ref. 20), Wei, Levinsky, and Su (ref. 18), and others. The complex per- 
turbation potential for this model is given by equation (A5). If the body can be ignored 
and if the vortex-sheet integration is performed in the physical plane, this equation can 
be written as 
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I 


where s is the distance along the wing and vortex sheet from the wing center line, and 
where the strength of a vortex sheet or wing segment is 

^ ) . . d Afcftsj (C2) 

ds 

Let a contour be drawn around the system composed of the wing, vortex sheets, 
feeding sheets, and vortex cores as shown in figure 10. With Cauchy's integral theorem 
the complex potential W can be written as 


w 


-C7” 


r(x,s) 

. X(s) - X 


ds 


a -*2 r(x,s) 


r s=j 
4s- 


y 2 


X(s) - X 


ds 


n 

s=y 


s ~ s max + 

2 x < s ) - x J S=-X V 


f S= ' x *( s max) AW(x,S) f S=X v AW(x,S) dg 

Js = -X.*. S - X ^S=X(s_„J s - X 


+ lim 

p-0 


1 w=oj 1 - 27 t w(x,X v + pe iw ) ie ia) 


nLv~ 


pw= w„-2?r w(x,-Xy + pe iaJ )ie ia) 

du, + p 2 V -- V - da) 

pe iw -X J w=w 2 -Xy + pe l<0 - X 


(C3) 


where AW is the difference of the complex potential on the inside and outside of the 
feeding sheet. (See fig. 10.) 


-X*(s m ax) ^ s max^ 



Figure 10.- Contour around wing, vortex sheets, feeding sheets, and vortex cores. 
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It can be shown that for practical purposes, only the first three terms in equa- 
tion (C3) need to be retained. Consider the last two terms in equation (C3). With equa- 
tion (Cl) the first of these terms can be written as 


lim 

p—0 


pa>=u>r27T w(x,X v + peH ie iw 

p \ — — : dco 

J U= X v + pe lU) - X 



(C4) 


where it is assumed that 


|x v - x| » p 

It can be shown that to the lowest order in p 

1 /( x v + pe iw ) - y| = ft - yi + 

|/ X v " ^2 


Consequently, equation (C4) can be evaluated as 


im p \ 

)— 0 ° w= u> 


lim 
P 


tOj-277 w(x,x v + pe ia5 )ie ia> 


X v + pe ia> - X 


du> 


a)= a). - 2tt 


2n 


(X v - X) 


n a>= O). 

lim p \ 1 

p—0 coj 


log. 


XyPe 




/X? - yh 


le* 40 da> = 0 


In a similar fashion it can be shown that the last term in equation (C3) vanishes. Consider 
the fourth and fifth terms in equation (C3). It can be shown that AW has the values r v 
and - r y on the feeding sheets from X^S max ) to X v and -X*(S max j to -X* , 
respectively. It follows that 


f s=x v AW(x,S) 

^S=x(s ma3 ^ S * X 


ds = -r v log. 


l - 



(C5) 
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I 


S=-x£ 


s=-x 


( s max) 


AW(x,S) 
S - X 


dS = T v log e 


; x;-x*(s max j 

A u. 


X* + X 


Assume that several turns of the vortex sheet are included as shown in figure 11. 
this figure it can be seen that for points X located near the wing surface 


X v - x(s max ) 


« 


x v - X 


Xy - X*(s max j 


« 


Xy + X 


X-plane 


Center line 


N 


JL 

/ / X v \ \ 

\\ \ '/ 
\ ^ - •> / 1 
S'/ 




(a) Physical plane. 



(b) Vortex- sheet strength. 

Figure 11.- Variation of vortex-sheet strength. 


(C6) 

From 
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It follows from these inequalities and equations (C5) and (C6) that 


C S=X V AW(jg,S) d „ 

J s-x(s max ) s - X 


Xy “ X(s max ) 


I 


S--X* AW(x,S) 

S - X 




dS 


Xy - X 

1 

x v - x | • 

- T — 

r 

x£ - x*(s max ) 

Ay - 

A V 

x£ + X 


. « 1 


. « 1 


It is concluded that the potential for separated leading-edge flow can be approxi- 

mated as 


0i 


Re 


s=-y. 


±r- 

iri J_ 


r(X,s) 

s=-s max X(s) " X 


ds + 


±c s=y 

^s=-y 


2 zM ds + 1 


y 2 


S - X 


171 


r s - s max y (x,s) d= 
i J s= $ 2 X(s) - X 


(C7) 


where it is assumed that the flow model includes several turns of the vortex sheet as 
shown in figure 11. 


Consider the derivative <p'^. From equation (C7) it can be shown that 



y'(x,s) 


s - X 


1 

S - X 



+ 


y(x,s ma x) 


X’(smax) 

X(s max ) - X 


X* ( s max) 
X*(s m ax) " X* 


X ( s max) 

X (s m ax) + X 


X ( s max) 
X(s m ax) + X 



(C8) 


In the derivation of this equation, the facts that for a symmetric flow field y(x, s) is anti- 
symmetric about the wing center line and that the quantity y'(x,s) is very small at all 
points on the vortex sheet have been used. The last condition, which arises because the 
vortex sheet can support no force, is obtained from equation (20) of reference 18. It is 
well known that the magnitude of the vortex- sheet strength diminishes as distance along 
the sheet from the wing tip increases; a schematic of the type of behavior which Smith 
(ref. 20) found for a delta wing is shown in figure 11. Consequently, the contribution of 
the nonintegral terms in equation (C8) can be made arbitrarily small by increasing the 
size of the vortex sheet sufficiently. As a result, equation (C8) is approximated in this 
appendix as 
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0i(x,X,X*) 



y’(x,s)| 


s - X 


s - X* 


ds 


(C9) 


It should be noted that equation (C9) is exact for attached leading- edge flow. 


Solution for 02a 

This derivation is the same as that of reference 25. The general form of the par- 
ticular solution is given by equation (43). The derivative 0^(x,X,X*) is given by equa- 
tion (C9). As stated previously, the Kutta condition applies at the wing leading edge so 
that to lowest order 



(CIO) 


Equation (CIO) is exact if the flow is attached at the leading edge. From equations (43), 
(C9), and (CIO), the particular solution for 02 a can be written as 


= - y + 


^(f ? /(*, s) ds 
32 tt ^ Ws=-y 2 


t_ ^2 / (x,t) 

■y 2 s " 1 


dt 


X * log e^#j + X log e(^ 


- 2 f S 72 y*(x,s) log e (s - X) ds f * 72 y’(x,t) log e (t - X*) dt 
J s=-y 2 J t=-y 2 


(Cll) 


From equation (Cll) it can be shown that the normal derivative of 0^ at the wing sur- 
face is 


902a(x,y,±°) _ y + 1 q 
9z ^ 4 ax 


y’(x,y) P.v. f y 2 y'(x,s)( — + log Is - y|)ds 

\ s ~ y / 


s= -y 


(C12) 


where the notation P.V. denotes the principal value of the integral. 

The potential 0 2a must satisfy a homogeneous Neumann boundary condition at the 
wing surface. As a result, the general form of the complementary solution of 0 2a is 
given by equation (47). With equation (C12) the potential 0 2a can be written as 


■ C _ y + l C y 2 _2_. 
^2a 1677 Jo--?? 

s ~ y 2 


'(x,s) P.V. r y'(^>t)(rz~g + l°Se I s ■ tlj dt 
^2 


logg |(s - X)(s - X*)J ds 


(Cl 3-) 
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At large radial distances r from the configuration, equation (C13) can be written as 


■»2a = l0g « ? 11 If 


s= ■ 


'^2 


y'(x,s) P.V. ^ 2 y'(x,t)^-^-^- + log e |t - s|jdt 
^2 


ds (C14) 


Solution for 02b 

The derivation of this solution is the same as that of reference 25. The general 
form of the particular solution is given by equation (44). The appropriate form of the 
potential 0^ to be used in evaluating this solution is given by equation (29). It follows 
from equations (44) and (29) that 0^ near the wing surface and the normal derivative 
of 02 b at the wing surface are, respectively. 


= 

and 



d<t>2 b (%>7 >±Q) 


81 


m ± (x,y) g"(x,y) - m^(x,y) 1 - g’(x,y) 


(C15) 


(C16) 


The potential 02k must satisfy a homogeneous Neumann boundary condition at the 
wing surface. Consequently, the general form of the complementary solution of 02k * s 



x f 

877 


s=y. 


8n J s=-y 2 L 


30 


P 

2b 


(x,s,+0) 

dZ 


30^ b (x,s,-O) 

dZ 




(Cl 7) 


When equation (C16) is substituted into equation (C17), it is found that the complementary 
solution is 


02b = “ ^“^2 J ~ S"(x,s) - m’(x,s) [l - g’(2,s)jjlog e [(s - X)(s - X*)j ds (C18) 

^2 

where the function m is given by equation (49). At large radial distances r from the 
configuration, equation (C18) can be written as 



—^2 log e ? f (m(x,s) g M (x,s) 

2tt\ 2 e J s=-y 2 l 


m 



ds 


(Cl 9) 
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Solution for <p 2c 

The general form of the particular solution for $ 2c i s given by equation (46). 

With this equation and equation (29) for the potential 0, near the wing surface, it can 

P X p 

be shown that 0 2c near the w ing surface and the normal derivative of 4 > 2c at the wing 

surface are, respectively, 


<P 


P 

2c 


(x,y,z) 



(C20) 


and 


84>f c (x,y,±0) 

sz 


\ m ± (x,y) 


(C21) 


The potential 0 2c must satisfy homogeneous Neumann boundary conditions at the 
wing surface. Consequently, the general form of the complementary solution of the poten- 
tial $> 2 C i s 



94>f c (x,s,+0) 

dZ 


9^f c (X,S,-0) 


dZ 


log e [(s - X)(s - X* j| ds 


(C22) 


With equation (C21), equation (C22) can be written as 



i r s= y 2 

AtiX 2 ^s=-y 2 


m(x,s 



s - X)(s 


X*) 


ds 


(C23) 


where the function m(x,y) is given by equation (49). At large radial distances r from 
the configuration, equation (C23) can be written as 



1 r s=y 2 

— log e r \ * m(S,s) ds 

2rr\ 2 J s=-y 2 


(C24) 
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THIRD- ORDER INNER POTENTIALS 


Consider the third- order terms in inner expansion (eq. (16)). When this expansion 
is substituted into equation (7) subject to the assumption that a, and v are related 
by equation (54), it is found that the potential $3 3 is governed by the two-dimensional 
Laplace equation in the cross-flow plane and that the potentials $3 3 and <p^ j are 
governed by the Poisson equations 


^2^3,2 = (y + 1) 


I 

0lfe ^6,1 + 02, 2 ) 


(Dl) 




0lfe ^ + 02,1 


3 r 9 r p 2 80 80 ) 


(D2) 


It is also found that the potential 4 > g is governed by a Poisson equation with the 
same degree of complexity as equation (7). For large values of the inner radial vari- 
able r, the equation for <£3 can be written as 


£2 


v 203 = (x + 1 )(0f02 j + • • • (D3) 

The particular solutions to equations (Dl), (D2), and (D3) for r » y 2 are 


<P 


P _ _ Y + 1 
3,2 " ‘ ~ 2 


i( F e F e) + ^ 




r log- r sin 0 


0 3,i = 2l (y + 1) f ’fe g a + g 2,i 


k 2 f " _ y + 1 

n n 


2a ^ L 


f’(F e F; 


‘r log- r sin 0 


y + 1 


4a 2 L 


f'(F e F4 


r log 2 r sin 0 
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It can be seen that if the complementary solutions $g g and 0g g are of the form 



the sum of the third- order terms in the inner expansion when written in terms of outer 
variables for large values of the inner radial variable r satisfies the equation 


sin^ a 


| log e jfc ^3,3 + log e b 0 3,2 + log e b *3,1 + *3 


= O 


sin^ a 


v 


log e t 


= Ol 


l0 Sexi; 


« y 2 - 


It is seen that this sum is not of lower order than The same relationship can be 
shown to hold for the sums of the fourth and higher order terms of the inner expansion. 
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